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LONG-TERM GOALS 

The long-tenn scientific objective of the Quantifying, Predicting, and Exploiting (QPE) 
Uncertainty Directed Research Initiative (DRI) is to improve the assessment of uncertainty in 
observations and predictions of sound propagation in littoral regions. 

OBJECTIVES 

The goal of this research is to understand and exploit the effects of the ocean state on acoustic 
propagation and detection. This work will contribute to that goal through regional ocean modeling 
and data assimilation. The modeling will include forecast and predictability studies to see the 
growth of uncertainty in time and space and the predictability of the propagation conditions on the 
shelf north of Taiwan from the ocean state. 

APPROACH 

The QPE DRI is a coordinated effort in which many types of measurements have been collected 
over the continental shelf to the north of Taiwan. The field results will be used to help characterize 
the rapidly varying physical environment and to study acoustic propagation and scattering in the 
region. 

The technical approach centers on using the Frechet derivative (adjoint) of a regional ocean model 
to explore the sensitivity of critical environmental features of the shelf region to the ocean in the 
region, to forcing, and to boundary conditions. In the longer term, acoustic remote sensing data, 
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together with data from direct measurements and satellite remote sensing, will be assimilated into a 
regional ocean model to estimate the evolving ocean state using the adjoint technique. 

Initially we configured a fine-resolution (1/12°) regional general circulation ocean model, using the 
MIT general circulation model (MITgcm), for a region encompassing Taiwan. This model can 
provide lateral boundary conditions to a smaller regional model to the northeast of Taiwan being 
run by Lermusiaux (MIT). Initial and lateral boundary conditions for our model are obtained from 
a global 1/12 Hybrid Coordinate Ocean Model (HYCOM) assimilation product. The choice of 
lateral boundary conditions is central to the issue of uncertainty. HYCOM is forced with daily 
Navy Operational Global Atmospheric Prediction System (NOGAPS) fluxes and uses the Navy 
Coupled Ocean Data Assimilation (NCODA) to assimilate data. First we needed to determine the 
fidelity of the output in the larger northwest Pacific region by considering the representation of 
those processes that will be felt at the boundaries of our regional model. These are primarily 
westward propagating mesoscale eddies. We also needed to consider the variability of the 
surrounding waters such as the annual cycle, mean flows, stratification and mixed layer depths. 

In preparation for the regional modeling/data assimilative effort we gathered both in-situ and 
remotely sensed data for the larger northwest Pacific region. These data include temperature and 
salinity vertical profiles from Argo floats, Expendable Bathythennograph (XBT) temperature data, 
velocity and temperature at 15m from surface drifting buoys, sea surface height from altimetry, 
and sea surface temperatures. The data was quality controlled and manipulated as appropriate for 
insertion in the data assimilation scheme. 

Once the regional model and data assimilation machinery was working, sensitivity studies were 
conducted prior to the FY09 experiment to aid in the design of the experiment and to understand 
the sensitivity of the models to the various data types and geometries. The intent is to be ready to 
combine the data obtained during the FY09 experiment with realistic ocean models when the 
observations become available in order to provide accurate estimates of the ocean state. 

WORK COMPLETED 

In the past year, the project focus has turned from evaluation of existing global model output to 
regional modeling of the Taiwan region. Two sets of regional MITgcm runs have been conducted, 
one in a large region that extends from 116°E to 140°E and from 17°N to 27°N, and the other in a 
domain one quarter the size, which extends from 116°E to 128.5°E and from 22°N to 27°N. The 
larger domain has been run at 1/12° and 1/24° degree, while the small domain has been run at 1/24° 
degree. Both setups used 50 layers in the vertical, including 5m spacing near the surface. 

The simulations were evaluated based on the realism of environmental features seen as keys for the 
observational program. The Kuroshio path was compared to float observations by Centurioni and 
Niiler, including both the Luzon Strait Loop Current and the Kuroshio surface water intrusions 
onto the shelf northwest of Taiwan. The similarity of the "Cold Dome" feature in the model to 
historical observations was also examined. 

The simulations in the large domain were run with a variety of initial conditions, boundary 
conditions, and forcing, in order to evaluate the model skill and its sensitivity to the controlling 
inputs ("controls"). Simulation periods of one to five years were tried, depending on the forcing 
product availability. 
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Initial and boundary conditions from both 1° ECCO and 1/12° HYCOM were used, and the 
higher-resolution HYCOM boundary conditions seemed to produce a better Kuroshio behavior 
north of Taiwan. A number of forcing fields were tried, including the NCEP/NCAR Reanalysis, 
NOGAPS, QSCAT winds, a combination of NCEP fluxes and NOGAPS winds, and similar 
combinations of NCEP or NOGAPS fluxes with QSCAT winds. No forcing product was clearly 
superior, but the NCEP fluxes with QSCAT winds seemed to have good perfonnance. 

Three topography estimates were tried: ET0P02, HYCOM bathymetry, and NCOR bathymetry 
supplied by the Taiwanese collaborators. The NCOR bathymetry had the most realism in the key 
experimental region NW of Taiwan, but the HYCOM topography gave the best Kuroshio behavior 
across the Luzon Strait. Simulations with NCOR topography showed several large intrusions of 
the Kuroshio into the Luzon strait, not consistent with observations, also calling into question the 
behavior of the Kuroshio along the east coast of Taiwan after exiting the Luzon strait. Experiments 
with doubled resolution and a variety of perturbed boundary conditions (in addition to the forcing 
products) did not change these results. Although some sensitivity analyses were carried out in the 
larger domain, the results were subject to concerns about model resolution and Kuroshio path. 

The smaller domain runs were adopted to focus more on the experimental region, and to eliminate 
the issue of the Luzon Strait Loop Current. The smaller domain allowed doubled horizontal 
resolution, also a key component of realism in the interaction with topography. To take advantage 
of recent data, the simulations were run for recent years: 2004 through 2008. Both NOGAPS and 
QSCAT winds were used with NCEP fluxes, and both wind products showed reasonable features. 

The forward runs showed realistic Kuroshio position, including intrusions onto the shelf, and had 
generally good cold dome structure. Cold features were formed both by interaction of the 
Kuroshio with the topography and by strong wind forcing on the shelf. The cold features were also 
influenced by the Taiwan Strait flow and its seasonal variability. The cold dome timescale was 
seen to be approximately ten days. The upwelling of high-salinity Kuroshio water was seen in the 
center of the cyclonic eddy forming the cold dome in some cases. 

Several typical events were chosen for sensitivity analysis, running the adjoint model backward in 
time from the location and time of a cold dome or intrusion event. The sensitivity maps showed 
upstream dependence of the cold dome and intrusions in line with expectations, and showing both 
growth and decay of sensitivity with time lag, depending on location. These results were delivered 
to the Chief Scientist for use in planning the experiment. 

The secondary focus of the project this past year has been the continued analyses of the 1/12° 
HYCOM data assimilative and forward simulations in the northwestern Pacific and preparation of 
results for publication. The entire duration of both simulations have been analyzed in terms of the 
predictability of intrusions of the Kuroshio into the East China Sea and the regional circulation 
near Taiwan. Specifically we have examined how the remote oceanic signal in the form of 
westward propagating mesoscale eddies interacts with the Kuroshio causing it to often first 
meander offshore and then intrude over the shelf north of Taiwan. These results complement the 
downscaled findings obtained from the regional models that use HYCOM output as lateral forcing. 
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RESULTS 


The adjoint model has worked successfully to show the sensitivity of key environmental features 
on the upstream ocean state, information used for the design of the observational programs. The 
sensitivities will allow a rapid propagation of uncertainties forward in time to regions of interest. 

IMPACT/APPLICATIONS 

This study will lead to the improvement of the assessment of uncertainty in observations and 
predictions of sound propagation in littoral regions. 

TRANSITIONS 

Methodology and data results can be made available to Navy scientists. 

RELATED PROJECTS 

The work described here is in collaboration with Professor Peter Niiler at SIO and Dr. Pierre 
Lermusiaux at MIT. 
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